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We calculate the magneto-resistivities of a two-dimensional electron gas subjected to a lateral 
superlattice (LSL) of rectangular symmetry within the guiding-center picture, which approximates 
the classical electron motion as a rapid cyclotron motion around a slowly drifting guiding center. We 
explicitly evaluate the velocity auto-correlation function along the trajectories of the guiding centers, 
which are equipotentials of a magnetic-field dependent effective LSL potential. The existence of 
closed equipotentials may lead to a suppression of the commensurability oscillations, if the mean free 
path and the LSL modulation potential are large enough. We present numerical and analytical results 
for this suppression, which allow, in contrast to previous quantum arguments, a classical explanation 
of similar suppression effects observed experimentally on square-symmetric LSL. Furthermore, for 
rectangular LSLs of lower symmetry they lead us to predict a strongly anisotropic resistance tensor, 
with high- and low-resistance directions which can be interchanged by tuning the externally applied 
magnetic field. 
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I. INTRODUCTION 

Pronounced commensurability oscillations of the mag- 
netoresistance of a two-dimensional electron gas (2D EG) 
subjected to a perpendicular magnetic field and a lateral 
superlattice, now also known as Weiss oscillations (WO), 
have first been observed on systems with a periodic mod- 
ulation in one direction (lD).ETu Subsequent work on sys- 
tems with a 2D lateral superlattice (LSL) showed that the 
modulation in the second direction tends to suppress the 
commensurability oscillations .QH The WO observed on 
samples with a ID LSL, and their suppression in sam- 
ples with a 2D LSL, were first explained by quantum 
mechanically. A ID modulation broadens the Landau 
levels into bands of oscillatory width, with finite group 
velocity. This leads, in addition to the scattering induced 
magnetoconductivity (the so called "scattering conduc- 
tivity"), to a "band conductivity" which vanishes if the 
Landau bands become flat .Bum A modulation in the sec- 
ond lateral direction splits these Landau bands into nar- 
row subbands with small group velocities, and, as a func- 
tion of the magnetic flux per unit cell of the 2D LSL, 
a self-similar energy spectrum ("Hofstadter's butterfly") 
results. EJEl It has been argued that this subband split- 
ting leads to a suppression of the band conductivity, if 
the modulation-induced width of the Landau bands is 
sufficiently large—and collision broadening effects are suf- 
ficiently weak.l'Eil Experiments on samples with a weak 
2D modulation and not too high mobility show indeed 
commensurability oscillations very similar to those ob- 
served in ID LSLs, which are suppressed (and changed 
in character) with increasing modulation strength and 
mobility.^ 

Soon after their discovery, Beenakkenl3 explained the 
most prominent of the WO in an electrostatically de- 
fined ID LSL classically, as resulting from an oscillatory 
E x B drift of the guiding centers (GCs) of cyclotron or- 
bits, where the GC velocity plays the role of the group 



velocity in the quantum treatment .0 The GC picture 
can be justified for weak modulations and .intermediate 
strengths of the applied magnetic field.oli3 It has been 
used to calculate the resistivity for different analytical 
forms of electrostatic, magnetpatatic and mixed modula- 
tions defining ID or 2D LSLsJilJ A clear classical picture 
for the suppression of the WO in 2D LSLs has, how- 
ever, not been developed for nearly a decade, although 
direct numerical evaluations of the diffusion tensor on 
the basis ofj-,ciassical ballistic models indicated such a 
suppression .QO ._. 

Recently Grant et alS3 emphasized that the GCs move 
along the equipotential lines of an effective potential, de- 
termined by an average of the modulation fields over 
unperturbed cyclotron orbitsO They argued that, in a 
2D LSL of square symmetry, these equipotentials are 
closed, and that therefore the GC velocity averages to 
zero, resulting in a suppression of the WO, and they 
confirmed this conjecture by ballistic model calculation. 
Recent experiments demonstrate also, that an asymmet- 
ric 2D modulation leads to much stronger commensur.a-| 
bility oscillations than a square-symmetric one does.c2l 
While Grant et alS3 employed the GC picture to make 
the suppression of WO in 2D LSL plausible, they did 
not really use it as basis of their calculation. Moreover, 
they considered a strong modulation, so that their re- 
sults are not directly comparable with previous predic- 
tions for weak modulationlij A consistent evaluation of 
the GC approach, which is known to yield a very simple 
and intuitive picture of the WO in ID LSLs, is so far not 
available for the case of 2D LSLs. The aim of the present 
work is to fill this gap. 

In Sect. H we discuss the GC approach and its limi- 
tations. Based on a numerical evaluation of the diffusion 



tensor, we present in Sect. Ill simple analytical results 



for the conductivities of square and rectangular LSLs 
defined by harmonic electric and magnetic modulation 
fields. These results depend only on one (square LSL) 
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respectively two (rectangular LSL) parameters, which 
are determined by the (seven) modulation model pa- 
rameters, the mean free path, and the average magnetic 
field. Illustrative examples are given, including a strongly 
anisotropic case in which the applied magnetic field can 
interchange the directions of high and low resistance. 
Mathematical details are given in two appendices. 



II. THE GUIDING CENTER PICTURE 

A. Heuristic definition 

We consider a 2D EG in the x-y plane subjected 
to a strong homogeneous, perpendicular magnetic field 
Bo = (0, 0, Bq) and a LSL defined by weak electric and 
magnetic modulation fields. The classical magnetocon- 
ductivity of this system can be calculatedEJ from the mo- 
tion of electrons at the Fermi energy Ep = (m/2)vp. 
Within the GC picture, this is assumed to be the super- 
position 



r(t) 



(i) 



of a rapid cyclotron motion r cyc (t) = i?(sina, — cos a) 
around a slowly moving guiding center r gc (i), where 
a(t) = u>ot + ao describes a uniform circular motion 
with cyclotron frequency luq — eBa/(mc) and radius 
R = vf/luq. This is obviously correct in the absence 
of modulation fields, where the position of the GC r gc (t) 
is a constant of motion, and in the presence of a homoge- 
neous in-plane electric field E _L Bo, where the GC moves 
with the constant drift velocity r gc = c(E x B )/Bq. 
For a perturbation by a position-dependent in-plane elec- 
tric field E = W(r)/e or perpendicular magnetic field 
B m = (0, 0, B m (r)), the GC picture is only approxi- 
mately valid, and several definitions of a "guiding cen- 
ter" are possible, which become equivalent in the limit of 
small perturbations. 

A reasonable candidate is the center of the circle of 
curvature at the point r(t). Taking the energy conser- 
vation (m/2)v 2 + V(r) — Ep into account and writ- 
ing the velocity as f = v = v(r) (cos ip, sin tp, 0\— with 
v(r) = vf[1 — V(y)I 'Ep} 1 / 2 , this center is given byc3 



r M = r + e z x v/(w + w mod ) , 



(2) 



where w mod = uj m + [e z x (v/w)] ■ Vu(r) with uo m (r) = 
eB m (r)/(mc). 

To lowest order in the modulation strength one may 
neglect the modulation effect w mo d in the denominator of 
Eq. (||). Then Newton's equation mv e[E + (v/c) x 
(Bo + B m )] yields the time derivativesEJ3 

E y u> m . E x Lo m , . 
xm = c— v x , y M = -c— v y . (3) 

In the spirit of the GC picture we may average Eq. (||) 
over the rapid cyclotron motion, i.e., we approximate v sw 



f cyc , replace r(t) = r gc (t) + r cyc (t) in the arguments of 
E and u> m , and take the average with respect to a = 
wot over one period. We assume that V(r) and w m (r) 
are periodic, with vanishing average values, on the same 
rectangular lattice with lattice constants a x = 2~k/K x 
and a y = 2ir/K y , 



V(t) 



q#0 



w m (r) 



to a e 



q#0 



(4) 



where q = (n x K x ,n y K y ). The averages over the cy- 
clotron motion can_be performed for each Fourier com- 
ponent separately.Ej The result of this approximation is 
the equation of motion for the GC, 



f 



'gc 



■ e z x Vy cff (r gc ) , 



(5) 



where the effective potential V e s (r) = X) q ^o e ' q ^q ff j * s 
determined 



= V^qR) 



mvF 



LOqJl(qR) 



(6) 



with Bessel functions Jo an( l Ji- According to Eq. (pj), 
the GC moves along the equipotential lines of the effec- 
tive potential T4ff(r gc ). Note that, in Eq. (H), we have 
identified the GC with the average value (over one cy- 
clotron cycle) of rjvf(i) as defined in Eq. (g), and that 
we have taken into account only terms in lowest order of 

Wmod 



B. Examples and limitations 

We have, for a large number of examples, integrated 
Newton's equation numerically to obtain the exact tra- 
jectories r(t) and rnj(t) as defined in Eq. (g), and also 
with the approximation ~^ 0. We found that 

in all cases with sufficiently small modulations (roughly 
l^modl ^ 0.2 wo) the average of r&f(i) and of its approxi- 
mation for w mo d — > are practically identical. Moreover, 
these averages agree with the average of the exact tra- 
jectory r(t) over the cyclotron motion, which we have 
calculated as 



r(i) 



T+(t)-T-(t)J T _ (t) 



T+(t) 



dt' r(i') 



(7) 



where ip(T±) = y(t)±7r, and tp(t) defines the direction of 
the velocity at time t, r(t) = u(r(i))(cos <p(t), sin <p(t), 0). 
This rather complicated definition of the time average 
seems necessary, since in the 2D LSL the velocity vec- 
tor f (t) is not an exactly periodic function of time, in 
contrast to the case of drifting orbits in a ID LSL.Iij 

As a typical example we show in Fig. [TJ, for two dif- 
ferent modulation strengths of a square-symmetric elec- 
trical modulation, rosette- like orbits together with rjifW 
(in the limit w mod — > 0) and r(i) defined by Eqs. (0) 
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and (Q), respectively. For weak modulation, both defi- 
nitions yield trajectories close to equipotentials, as ex- 
pected from Eq. (||). However, rjvf (i) exhibits rapid fluc- 
tuations around its equipotential, with an amplitude that 
increases with the modulation strength. Using the cy- 
clotron motion as a reference, we see from Fig. [| that 
the velocity of the GC motion increases (essentially lin- 
early) with increasing modulation strength. 
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FIG. 1. Parts of rosette-like orbits (about 39 cyclotron cy- 
cles, same initial conditions) in the electric modulation po- 
tential V(x,y) — e Ef[cos qx + cosqy], with q = 2n/a and 
qR = 2, for e = 0.05 (left) and e = 0.15 (right). Thick solid 
lines show the corresponding r(t) as defined by Eq. ([?]), thin 
dotted lines the equipotentials of V(x,y). The corresponding 
rA/(i), Eq. (Q) with u mo d = 0, are shifted downwards by one 
lattice period. 



acteristically from the equipotential trajectory predicted 
byEq. (§). 

The failure of Eq. (|) near the "flat-band conditions" 
is not very important for the calculation of the conduc- 
tivity, since there the GC drift is anyway slow and thus 
contributes little to the conductivity. There are, however, 
natural limitations of the GC picture. Of course, if the 
(average) magnetic field Bo becomes too strong, the clas- 
sical approach fails and Landau quantization effects must 
be taken into account. If Bq is zero or very small, "chap* 
neled orbits" occur similar to the case of a ID LSL.E3 



In contrast tc 
chaotic orbitsfc 



ID case, in 2D LSL one also observes 
1 21 l if the modulation is sufficiently strong 



and B is sufficiently small, w mo d > ujq. In the regime of 
chaotic orbits (i.e., for the model of Fig. [I] with e ~ 0.05 
for qR > 7) the GC picture is not useful, since it does 
not simplify the description of the electron motion. 
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If the modulation has only a rectangular instead of 
a square symmetry, the GCs may follow either closed 
(localized) or open (drifting) equipotentials. A typical 
example with an electrostatically defined LSL is shown 
in Fig. |. 

We want to mention that there are situations in which 
the GC picture works, but Eq. does not. This is, e.g., 
the case, if due to the Bessel functions in Eq. (^) the effec- 
tive potential vanishes for finite modulation, and there- 
fore the first order approximation Eq. (|^) fails. For the 
square-symmetric harmonic electric modulation consid- 
ered in Fig. this happens for Jo(qR) = ("electric flat- 
band condition"). Then, besides orbits with GCs mov- 
ing around potential maxima [r = (ma, na)] and minima 
[r = (2m + l,2n+ l)a/2], there are also orbits with GCs 
moving around saddle points at [r = (2m, 2n + l)a/2] 
and [r = (2m + l,2n)a/2], which are not described by 
Eq. (g). The approximation Eq. (||) becomes also poorer 
with increasing modulation strength. Thus, we see from 
the right thick line in Fig. [I] that the GC deviates char- 



FIG. 2. Thirty-eight cyclotron cycles of a drifting (left) and 
of a localized (right) orbit in the electric modulation poten- 
tial V(x,y) — 0.1 _Ef[cos ga; + 0.25 cosqy] (equipotentials indi- 
cated by thin dotted lines), with q = 2-K/a and qR = 2. Thick 
solid lines: GC trajectories as defined by Eq. (|). The rapidly 
fluctuating trajectories of rjf(t), Eq. (El) with cj mo d = 0, are 
shifted downwards by one period. 

In the following we will consider the decomposition of 
Eq. (Q) for the electron motion in a LSL and calculate the 
GC motion from Eq. (|). This is a good approximation 
if the modulation fields defining the LSL are sufficiently 
weak and if the average magnetic field Bq is sufficiently 
strong. We want to emphasize that the GC picture yields 
reasonable results for the magnetoconductivity even in 
situations in which Eq. ([!]) with r gc (i) calculated from 
Eq. (j^) does not yield a reasonable approximation for 
an individual trajectory r(t) with the same initial condi- 
tions. 
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C. Conductivity from GC motion 

To calculate the magnetoconductivity of a 2D EG in 
a LSL within the relaxation time approximation, we use 
Einstein's relation er M „ — D fJiU e 2 m/(iTh' 2 ) and the Cham- 
bers formulac3il3 for the diffusion tensor , which con- 
tains the velocity autocorrelation integral along a tra- 
jectory, averaged over all initial conditions r(0) = ro, 
f(0) = v(r ) (cos</?o,sin<£ )- With the decomposition 
of Eq. (0) this yields three types of terms. One term, 
which contains only the cyclotron velocity and must be 
averaged over the initial value ao, yields the Drude con- 
ductivity tensor. The mixed terms, containing both the 
cyclotron and the GC velocity, vanish upon averaging 
over ao- Finally, the term containing only the GC drift 
contribution is given by 



dte-^(vf(t)vf; c (0) 



init 



where the average has to be taken over all possible initial 
positions r gc (0) of GC trajectories in a unit cell of the 
periodic potential. 

With the dimensionless coordinates £ = K x x and rj = 
K y y and the effective potential 



(mK x ,nK v ) 



/Vcha j 



(9) 



£(t) = £(0), and v y — w'(^(0))/K y , independent of time. 
Thus we can immediately evaluate Eq. (0) to obtain 



yy 



D sc 

xy 



siv r d$ 



K 2 



.2tt 



and 



K x V c ha 



muio 



w'(0 



(12) 



For simple harmonic modulations, e.g., Vq = Sq^±K.o)VK 
and w q = 8q.(±K.o) UJ ±K with uj-k = u* K in Eq. (|J), one 
obtains from Eq. (^) V c s (r) = Vo cos(Kx + a), with Vq = 
2\V?3^\- With y cha = V and «/(£) = -sin(£ + a), 



Eq. ( |l2| ) reproduces the known formulalij 



1D e in as / 
yv = 1^2K 2 ' 



(13) 



(8) where the WO result from the oscillatory behavior of 



SI 



K 2 v 2 F 

UJ 



(14) 



Here a complex ratio lvk/Vk allows to describe a phase- 
shift between electric and magnetic modulation. 



B. Weak 2D modulation, Sir < 1 



where V c ha is an energy characteristic for V c s (r), e.g. its 
maximum, Eq. (|^) reads 



d£ dw d-q dw 
dt = drj' dt = ~ ~dZ' 



with 



SI = K x K y 



Vc 



oka 



mujQ 



(10) 



(11) 



Given the analytical form of the effective potential 
w (£,iV)> we can use Eqs. ( |lO| ) to calculate the diffusion 
tensor (^) as an integral over the equipotentials of u>(£, 77) 
(see appendix A). The results then depend only on the 
parameter Sir. In the following section we will demon- 
strate this with a few explicit examples. 



III. RESULTS AND EXAMPLES 

To keep the notation simple, we will in the following 
consider only superlattices with a rectangular symmetry, 
with ID and square-symmetric LSL as limiting cases. 

A. One-dimensional modulation 

If the periodic potential depends only on one coor- 
dinate, say w(£), Eq. (|io|) yields v x = £/K x = 0, 



For a 2D superlattice potential one obtains a similar 
simple result, if the modulation amplitude (or the relax- 
ation time r) is sufficiently small, so that Sir <C 1. Then 
we may approximate v^(t) ^(0) in Eq. (|^), so that 
the t integral becomes trivial, with the result 



d£ / dr\ W/i w p 



(15) 



with <t x = 1, a y = — 1, w x = dw/d^, w y — dw/drj and 
the notation ft = y (or x) if fi = x (or y). The factor 
(2tt) 2 is the area of the dimensionless unit cell. 

If the periodic potential is additive, u>(£, rf) = toW (£) + 
w^(rj), the off-diagonal elements vanish, D^ y — = 
0, and the diagonal elements agree with those of the cor- 
responding ID modulations, Eq. (|l2]). 

This weak- modulation limit Sir < 1, in which the 
magnetoconductivity is independent of the pature of the 
GC trajectories, has been discussed in Ref.li3. However, 
with increasing modulation strength (and larger relax- 
ation time, i.e. larger mean free path) the nature of the 
GC trajectories will become important. For Sir ^> 1 the 
time integral will be proportional to the average velocity 
along the trajectory. For closed trajectories, this average 
will vanish, whereas for open equipotentials, which may 
exist either in x or in y direction, the average may be 
finite. Thus we expect that, in the limit Sir — > 00, closed 
equipotentials will not contribute to the diffusion tensor, 
whereas the contribution of open ones will be similar to 
the case of ID modulation. 
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C. Square-symmetric harmonic modulation 

We now consider the 2D version of the simple har- 
monic modulation discussed in Sect. Ill A, i.e., assume 



in Eq. 



V a 



(<* q ,(±ic,o) + tiq,(o,±K)WK and uj. 



(<*q,(±K,0) + <*q,(0,±K)) w ±K witn U-K 

effective potential has the form 



to 



K ' 



Then the 



V eS (r) = V [cos(Kx 



cos(Ky + a)] 



(16) 



with an (irrelevant) phase shift a, and all equipotentials 
are closed lines around cither a maximum or a minimum, 
except those for V e g(r) — 0, which are straight lines. We 
find that the angular velocity of the GC drift along the 
equipotentials is given only by the parameter Q defined in 
Eq. ([Til), and geometrical factors. As a consequence, the 
suppression of the GC contribution to the conductivity 
can be expressed by a function $(f2r), and instead of 
Eq. (H) we obtain 



2D 



2D 

mi 



nh 2 2K* 



(17) 



and a*y = <7~~ = 0. The actual calculation of $(f2r) is 
sketched in Appendix 0. The numerical results are plot- 
ted as diamonds in Fig7]|, together with some analytical 
approximations, which are obtained from the asymptotic 
behavior of the correct result for small and large values of 
VLt (see Appendix |b|). Apparently the three-parameter 
interpolation formula 



$ 3 (ftT) = [l + 0.25(ftT) 2 ]/[l + 0.75(fh-) 



-0.076 (fir) 4 ] 
(18) 



provides a very good fit to the correct numerical result 
for all values of Qt. Note that $(Or) — * 1 for VLt — ► 0, 
as we expect for the weak-modulation limit. 



^(Hx), numerical 

* w cnx)=i/[i+(otr 

* 1 (nt)=l/[l+a(£2x)i 
*,(£2t) 




FIG. 3. Suppression of conductivity in 2D superlattice with 
square symmetry; numerical result $(f2r) (open diamonds), 
quadratic approximation with cut-off energy e qu = 0.228 
(long-dashed), one-parameter interpolation with a = 0.304 
(dashed), and the three-parameter interpolation of Eq. ( |li| ) 
(solid line). 



We want to emphasize that, for the square-symmetric 
harmonic cosine modulation the suppression of the GC 
induced contribution to the conductivity is described by 
the single parameter fir which, according to Eq. (M), 
itself depends on modulation strength and period, and 
on the cyclotron radius R — vf/loq. 
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FIG. 4. (a) GC drift contribution to the conductivity 
versus magnetic field in units 1/qR for electric modulation 
V(x,y) — 0.02i5.F[cos(p; + cos qy] and several values of the 
mean free path between qX = 50 (lower dash-dotted) and 
1000 (upper solid line), (b) Same resistance data divided by 
(qA) 2 , compared with result for the one-dimensional modula- 
tion V(x) = 0.02Ef cos qx (thick line) which is independent 
of q\. 

As an instructive example, we plot in Fig. |](a), un- 
der the assumption luqt 1, the GC contribution 
Apxx/po ~ {^qt) 2 &<Tyy / <tq for the electric modulation 
V(x, y) — 0.02_Ep[cos qx+cos qy] and several values of the 
mean free path A = VpT (ao = 1/po = e 2 n c \T/m). Away 
from the flat band conditions, given here by the zeroes of 
the Bessel function Jo, fir becomes large with large mean 
free path, and Ap xx /po w (Q,t) 2 Q(Q,t) / (qR) 2 approaches 
the limiting curve 3.29 /(qR) 2 , which is also indicated in 
Fig. |(a). We see that, as compared with the Drude re- 
sistance po of the homogeneous system, the modulation 
induced correction to the resistance increases with in- 
creasing mean free path and finally saturates. Since the 
shape of the resistivity curves in Fig. [| depends only on 
the parameter fir and is proportional to the modula- 
tion strength, variation of the modulation strength leads 
to a set of curves similar to that shown in Fig. |](a) for the 
variation of (qX) 2 . With increasing modulation strength 
the curves will saturate and approach the same limiting 
curve, indicated as thick dash-dotted line in Fig. ||(a). 

On the other hand, if we compare this modulation cor- 
rection calculated for the square symmetric case with 
that obtained for the corresponding ID modulation, we 
find, with increasing mean free path, an increasingly 
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strong suppression of the WO. This becomes evident from 
Fig. [f|(b), where we have divided Ap xx /po by (</A) 2 , since 
with this normalization the ID result becomes indepen- 
dent of qX. In this plot the suppression of the resistivity 
maxima becomes stronger with increasing mean free path 
and increasing modulation strength. 




1/qR 

FIG. 5. Comparision of resistivities calculated for electric 
modulation V(x,y) = 0.02Ef [cos qx + cos qy] within the GC 
picture (solid lines) and from numerical solution of the lin- 
earized Boltzmann equation (dotted lines), for two values of 
the mean free path. 

To show that the GC picture yields reasonable results, 
we compare in Fig. || calculations based on Eqs. jl7| ) 
and ( |is| ) with results obtained numerically from the lin- 
earized Boltzmann equation. The Boltzmann equation 
was solved by Fourier expansion of the distribution func- 
tion with respect to the periodic position variables x and 
y and the angle tp in velocity space, similar to the pro- 
cedure described in RefEj for a ID LSL. To obtain the 
curve for qX — 1000 with sufficient accuracy, about 40 000 
Fourier coefficients had to be included, and the calcu- 
lation, using an optimized parallel code, took about 6 
hours on a CRAY-T3E with 128 nodes. The compari- 
son shows that the GC approach with the approximation 
(|l8|), which requires only negligible numerics, yields sur- 
prisingly good results for weak modulations. The agree- 
ment will become poorer for stronger modulation and for 
much smaller qX. Then, with decreasing £?o, the max- 
ima of the Ap xx oscillations in the GC approach will 
still extrapolate to zero, whereas the correct calculation 
yields damped oscillations around a nearly constant, fi- 
nite Ap xx value. But this difference occurs also for ID 
LSLs and is well understoodJlj 

We conclude that the GC approach yields reasonable 
results for not too strong modulations (and not too small 
qX values), and we will use it as a versatile approach to 
discuss interesting situations of lower symmetry. 



Pure magnetic modulations lead to similar results as 
pure electric modulations, of course with modifications 
due to the differences between the Bessel functions J\ and 
Jo, notably a phase shift. Interesting new situations oc- 
cur for mixed electric and magnetic modulations, which 
can be achieved experimentally, e.g., by bringing a rect- 
angular pattepi of ferromagnetic islands on the surface 
of the sample.Ej Superpositions of harmonic electric and 
magnetic modulations, eventually with a phase shift, can 
easily be evaluated using Eq. ([18]) , provided the effective 
potential according to Eq. (|6|) has square symmetry. 



D. Harmonic LSL with rectangular symmetry 

A LSL with exact square-symmetry is an idealized lim- 
iting case and hard to realize experimentally. Therefore, 
we now consider the more general case of a rectangular 
LSL, which allows to interpolate between ID and square- 
symmetric 2D modulations, and to approach both limit- 
ing cases. To keep the discussion simple, we restrict it 
on harmonic electric and magnetic modulations in both 
directions, so that the effective potential is of the form 

V eS (r) = V x cos(K x x + a x ) 4- % cos(K y y + a v ) , (19) 

where the ratios of amplitudes and phases may depend on 
the amplitudes and relative phases between the electric 
and magnetic modulations in x and in y direction, and, 
in contrast to Eq. (0), on the average magnetic field Bq. 

Besides its simplicity, this model is important for the 
physical reason, that higher modulation harmonics de- 
crease exponentially with the distance of the 2D EG from 
the surface if the modulation is produced by some type of 
surface structuring. Thus, if this distance is large enough, 
it will be sufficient to consider only the basic cosine mod- 
ulation. 



1. Numerical and analytical results 

For a given modulation, the ratio V y /V x in Eq. ( |l9| ) 
may change magnitude and sign as a function of Bq . This 
can lead to interesting switching effects which we will 
discuss below. For the calculation of the conductivity 
components (Appendix ^), we assume however always 
< k — V y /V x < 1, which may eventually require an 
interchanging of x and y in the final results. Then, with 
V c hn = V x in Eqs. @ and (pi]), and with a suitable choice 
of the origin, the dimensionless potential (S) becomes 



w (^v) = cos £ + K cos V i < k < 1 . 



(20) 



For k — we have the ID modulation in x direction, 
and the equipotentials are straight lines in y direction. 
For k = 1 we have the square-symmetric case where all 
equipotentials are closed lines. These cases have been 
considered above. For < K < 1, there exist closed 



G 



equipotentials with 77) = e around maxima in the 
energy interval 1 — n < e < I + n, closed equipotentials 
around minima in the interval — (1 + ft) < e < — (1 — «), 
and open equipotentials in y direction for — (1 — ft) < 
e < 1 — k. We can show that, in the limit of large mean 
free path (r — > 00), the GC contribution Aa xx comes 
only from closed orbits, and shows a suppression similar 
to that obtained in the square-symmetric case. The con- 
tributions to Aa yyi on the other hand, come from both 
closed and open equipotentials. The latter lead to an 
increase with increasing r, similar to the ID case. 



The off-diagonal components Aa xy = Aer 



can 

be shown to vanish from symmetry reasons. The analyt- 
ical considerations of Appendix [A] show that the diagonal 
components can be written as 



A(7 W = 



(21) 



limit $ aa (cx), ft) is easily calculated numerically and well 
approximated by 



$ 2 ( K ) = 1 - 1.645 k 3/2 + 0.645 n 5/2 , 



(23) 



see Appendix [B|. Incorporating this into an interpolation 
formula that reduces for ft = 1 to the previous fit (pT 
we obtained 



$( 3 )(Slr,«) = $ 2 ( K ) + 



[1-$ 2 (*Q][1 + c*k(^) 2 ] 
1 + K+/3 k )(S1t) 2 + 7k (^) 4 ' 

(24) 



with a K = 0.25 sin 2 (tw/2), (3 k = 0.5k 2 /[1 - $ 2 («)], and 
7 K = 0.076 sin 2 (7TK/2). This approximation is indicated 
by the lines in Fig. |(] and will in the following be used 
instead of $ TO (fir, ft). 



with ft — K x K y V x /(m(jj ), so that the (suppression) ef- 
fect of the 2D modulation now is described by two pa- 
rameters, Sir and ft. We have numerically calculated the 
functions $ w (J1t, ft), which of course satisfy the consis- 
tency relations Q xx (£It, 1) = $ yy (riT, 1) = $(f2T). Since 
it is rather time-consuming to calculate the successive 
fourfold integrals with sufficient acccuracy for each spe- 
cific example anew, we tried to fit the <&p M (ilT, k) by sim- 
ple analytic expressions. We found that the numerical 
results for <j> xx (jlr, k) are very well (with an error of less 
than 1 per cent) approximated by k 2 $(S1t), so that a 
good approximation is 



1.0 4 



o 0.4 




FIG. 6. Numerical result for 5> TO (fir, ft) (diamonds) for 
several values of the anisotropy parameter k = V y /V x . The 
solid lines are for the approximation <& y 3 y , Eq. (p^). 

$ ra (Slr, ft) « k 2 $(S1t) » k 2 $ 3 (S1t) , (22) 

with $ 3 (ilr) defined by Eq. (|l^). Numerical results for 
$> yy (p,T, ft) are shown as diamonds in Fig. |^. Apparently, 
for ft < 1 they approach a finite limit for Sir — ► 00. This 



2. Tuio examples 

First we consider in Fig. [7] a purely electrostatic mod- 
ulation on a square lattice, a x = a. y = 2ir/q, V(x 7 y) = 
V x cos(qx) + V y cos(qy), so that = V l _ l \Jo(qR)\ and the 
ratio ft = Vy/Vx — V v /V x is independent of Bq, and 
n = q 2 V x \J (qR)\/{mu; ). 




% 1.5 



0.1 0.2 



0.3 0.4 0.5 0.6 

1/qR 



FIG. 7. GC contribution to the resistivities for the elec- 
tric modulation V{x,y) / Ef = 0.02 [cos(g:r) -\-ncos(qy) \ and 
qX = 400 . For k — 1 (thick line) Ap XI = Ap yy . For k < 1, 
the result for Ap yy lies below, that for Ap xx above this thick 
line, and both are indicated by the same line style. 

For k<1, there exist open equipotentials only in y di- 
rection. With decreasing k their number increases, and 
Ap xx increases towards the results for the ID modulation 
(ft = 0). Simultaneously Ap yy cx k 2 decreases, and van- 
ishes in the ID limit. The degree of anisotropy increases 
with both the modulation amplitude and the mean free 
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path, since, for fir ^> 1, Ap yy /p sa 3.29k 2 /{qR) 2 sat- 
urates, while Ap xx /p ~ (f2r) 2 $ 2 (K)/(qi?) 2 increases 
without limit. 

The anisotropy parameter n is only a constant inde- 
pendent of Bo if we have either a pure electric or a pure 
magnetic modulation on a square lattice, i.e., with the 
same period in x and y direction. In all other situations, 
the Bessel functions in Eq. (||) lead to a Bo-dependent k. 
In such cases we use the following convention to express 
the relevant parameters Q and k in terms of the original 
parameters specifying the modulation. 

We measure energies in units of Ep = mvp/2 and 
the average magnetic field in dimensionless units l/(qR), 
where R = vf/ujo is the cyclotron radius and q = 
yj K x K y . Then, for a suitable choice of the coordinate 
system, the modulation may depend on the following 
seven parameters: (1) the ratio of the lattice constants 
o. y /a x = K x /K y , (2) the amplitudes e v — V v /Ep of 
the electric cosine potential V(x,y) — Ep[e x cos(K x x) + 
e y cos(K y y)], (3) the amplitudes [i, v = 2uj v /(K v vp) and 
(4) the relative phases a v of the effective magnetic modu- 
lation potential Ep[p, x cos(K x x + a x )+p, y cos(K y y+a y )]. 
For each value of the average magnetic field, we can cal- 
culate from these seven model parameters the two pa- 
rameters of the effective potential (^|) which are relevant 
for the conductivity, namely the absolute values of the 
complex numbers e,jJo{K^R) + p v J\(K v R) exp(ia„), 

e v = ^\£ y J {K y R) + p v Jx{K v R) cosav] 2 

+ [p„J 1 (K„R)sma l/ } 2 y /2 (25) 

for v = x, y. The phases of these complex numbers can 
be compensated by a suitable shift of the coordinate sys- 
tem and have no effect on the conductivity. In the fol- 
lowing we use these two parameters in the form e max — 
man[£ x , e y ] and k = minle^ , eJ / max[e g , e v ] . Taking the 
characteristic energy in Eq. ( |ll| ) as V^ha = Emax-Ef, we 
obtain ft — a;o£ m ax(<7-R) 2 /2. 

To characterize the system completely, we have to 
specify the mean free path A = vpr, which we write in 
the dimensionless form qX, so that uiqt = X/R. Finally 
we obtain for the GC drift contribution to the conduc- 
tivity tensor 

= 71F2 I*-) 2 *W (^ 2 Ai?£max, , (26) 

with x = y and y = x if e m ax = £y, and with x = x 
and y = y if £ max = S x . Since in the regime of com- 
mensurability oscillations luqt 3> 1, the GC correction to 
the resistivity tensor is Ap^ lfl / po — (luqt) 2 Acr^/cro, with 
x, — y and y = x. 

As a very interesting example we consider a purely 
electrostatic modulation, but now on an rectangular su- 
perlattice with equal modulation amplitudes e x — e y but 
different periods in both directions, a y /a x — The 
interesting aspect of this model is that now the effective 



potential changes its symmetry as a function of the mag- 
netic field strength, since the arguments of the Bessel 
functions in Eq. (^|) are different. If one of the Bessel 
functions vanishes, i.e. if the "flat-band" condition for 
this direction is satisfied, the effective potential shows 
a purely ID modulation in the other direction. When 
the effective modulation potential in x direction is larger 
than that in y direction, there exist open equipotentials 
in y but not in x direction, and vice versa. Typical re- 
sults for the resistivity corrections are shown in Fig. ||. 
For relatively small mean free path as in Fig. |^(a), the os- 
cillations of the resistivity components Ap^^ look similar 
to the results one would expect for the corresponding ID 
modulations. At relatively low magnetic fields, there oc- 
curs however a kind of beating effect, manifested in a non- 
monotonous decrease of the oscillation amplitude of Ap xx 
(solid line) with decreasing magnetic field Bo. The reason 
for this non- monotonous Bo-dependence of the maxima is 
easily understood. The maxima occur nearly in the mid- 
dle between adjacent flat-band conditions Jo(K x R) = 0. 
If for these Bo-values the effective modulation in y direc- 
tion is large [i.e., if no zero of Jo(K y R) is close], the GC 
motion is essentially two-dimensional, and the maximum 
of Ap xx is suppressed below the corresponding one of a 
ID modulation in x direction. If, however, the maximum 
of Ap xx appears near a zero of Jo(K y R), the modula- 
tion in y direction is small, and the Ap xx maximum as- 
sumes a large value close to that of the corresponding ID 
modulation in x direction. This explains why the Ap xx 
maximum near (qR)^ 1 — 0.071 is higher than those near 
0.091 and 0.116. 
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FIG. 8. GC contribution to the resistivities for modulation 
V(x,y)/E F = 0.02 [cos(K x x)+cos(K y y)], with K x /K y = y/2, 
versus magnetic field in units of 1/qR, with q = */ K x K y ; (a) 
for q\ = 80, (b) for q\ = 800. The "fiat-band" conditions 
Jo(K x R) — and Jo(K y R) — are indicated by filled and 
open triangles, respectively. 

These anisotropy effects arc drastically enhanced for 
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larger mean free path, see Fig. ||(b). If, for example, at 
a given Bq value the effective modulation potential in 
y direction is smaller than that in x direction, the con- 
tribution to Ap yy (oc Aa xx ) results only from guiding 
centers moving along closed equipotentials, and Ap yy /po 
is bounded by 3.29(e y /e x ) 2 /(qR) 2 < 3.29/(gi?) 2 (cf. 
Fig. ^). There exist, however, open equipotentials in y 
direction which lead to an increase of Ap xx oc (qX) 2 with 
increasing mean free path. The result is a very effective 
switching as a function of the magnetic field Bo between 
regions with large Ap xx and very small Ap yy and regions 
with small Ap xx and large Ap yy , as is seen in Fig. ra(b). 

If one mixes electric and magnetic modulations with 
different phase shifts in both directions, one may achieve 
such switching effects also on a square lattice, K x = K y . 

IV. SUMMARY 

We have evaluated the modulation correction to the 
magneto-resistivity tensor of 2D EGs in LSLs of rectan- 
gular symmetry within the GC picture. We have em- 
phasized that this classical approach can be useful only 
within a restricted regime of sufficiently weak modula- 
tions and sufficiently strong (average) magnetic fields, 
where the electron motion may be approximated as 
a rapid cyclotron motion around slowly drifting GCs. 
Within this regime, in which a 2D EG with a ID LSL 
exhibits regular commensurability oscillations (WO), we 
have investigated the effects of the model parameters 
(modulation strengths, anisotropy, phase shifts) and the 
mean free path on the amplitudes of the WO. For har- 
monic electro- and magnetostatic modulations we have 
obtained essentially analytical results. 

The fact that the GCs move approximately along the 
equipotentials of a magnetic-field dependent effective po- 
tential, with a velocity essentially proportional to the 
strength of this potential, leads to an interesting depen- 
dence of the WO amplitudes on these model parame- 
ters. In contrast to a ID LSL, which has only extended 
straight-line equipotentials, a 2D LSL has also closed 
equipotentials around the extrema of the effective poten- 
tial. The difference between closed and extended equipo- 
tentials becomes important in the limit of a large mean 
free path, since the magnetoresistance is sensitive to the 
mean velocity of the GC motion between two successive 
scattering events. If the scattering time is sufficiently 
large, the GC velocity along closed equipotentials aver- 
ages to zero, and then these equipotentials do not con- 
tribute to the magnetoresistance, whereas the contribu- 
tion of extended equipotentials becomes very large. 

This leads to strongly anisotropic resistivities, if the 
effective potential has rectangular but not square sym- 
metry, and to very interesting magnetic-field dependent 
switching effects if the symmetry of the effective potential 
changes as a function of the average magnetic field. 

For the 2D EG with a weak square-symmetric mod- 



ulation we find with increasing mean free path an in- 
creasing suppression of the WO amplitudes below those 
obtained for the corresponding ID modulation. This re- 
sult provides a classical explanation of the suppression 
of the band-conductivity observed in early experiment^ 
on holographically modulated high-mobility samples,cl 
which previously has been explained with quantum argu- 
ments based on the-.subband splitting of the Hofstadter 
energy spectrumJjO For fixed mean free path, our (basi- 
cally analytical) result reduces in the iimit of very weak 
modulations to the predictions of Ref.Ej, without notice- 
able suppression of the WO. For realistic values of mod- 
ulation strength and mean free path, our present results 
yield, however, a strong suppression. We want to point 
out that our present classical explanation of the WO sup- 
pression and the previous quantum one are not contra- 
dictory. Both need for the explanation of an effective 
suppression a sufficiently strong modulation and large 
mean free path (i.e. weak disorder). 

Qualitatively our result is also in agreement with the 
recent prediction of the suppression of WO by Grant et 
al.t3, which applies to the case of intermediate mean free 
path and strong 2D modulation. From our investigation 
of trajectories we expect, however, that for this strong 
modulation the regime of small and intermediate values 
of the average magnetic field is dominated by chaotic 
motion, so that the GC picture cannot be expected to 
yield quantitatively correct results. 

Finally we want to comment on the fact that the 
Chambers formula (^J) contains a scattering time, which 
describes isotropic impurity scattering, whereas calcula- 
tions for the 2D EG with a ID LSL based on Boltzmann's 
equation have revealed that predominantly small-angle 
impurity scattering has to be considered for a quanti- 
tative understanding of the WO amplitudes. We did 
not try to go beyond the simple relaxation time ap- 
proximation in the GC picture, since (i) on the level of 
Boltzmann's equation, where we know how to describe 
anisotropic scattering, we cannot separate the GC from 
the cyclotron motion, and (ii) small-angle scattering of 
an electron between locally nearby trajectories may in- 
clude large changes of the corresponding GCs, and we do 
not want to introduce unjustified assumptions on scat- 
tering between GCs. In view of the general limitations of 
the GC picture, we rather want to consider the relaxation 
time r as a phenomenological parameter, which may be 
chosen to fit experiments qualitatively. We think, how- 
ever, that t should be considered as the total scattering 
time, which in the case of strongly anisotropic impurity 
scattering is much shorter than the transport or momen- 
tum relaxation time. 
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APPENDIX A: CORRELATION INTEGRALS 
ALONG EQUIPOTENTIALS 

To study the effect of closed equipotentials, we as- 
sume that the effective potential wl^rj) has either iso- 
lated maxima or isolated minima, or both. For instance, 
the model (|2^) has, for arbitrary integers m and n, iso- 
lated maxima at (£,77) = 2ix(m 1 n) and isolated minima 
a t (£,,v) = (2wi + l,2n + l)ir, and all equipotentials 
w (€) v) = e f° r |e| > 1 — « are closed. 

We assume that closed equipotentials around a max- 
imum (minimum), which we take as origin, exist in the 
energy interval e max > e > e sup (e min < e < e inf ). In 
terms of polar coordinates, 



£, = p cos ip, 7] = p sin ip, — 7r < tp < 7r . 



(Al) 



the equipotential with energy e is described by the equa- 
tion p = p e (ip), which maps ip onto the solution p of 
w(p cos ip, psmip) — e for fixed e and ip. Along the 



equipotential with energy e the equations (10) reduce to 



dip/dt = ±n/J e (ip) 



(A2) 



where the upper (lower) sign stands for orbits around a 
maximum (minimum), and 



COS ip + sin ip Wri 



(A3) 



We can use Eq. (|AJ) to substitute in Eq. (||) the inte- 
gration variable t by tp. Writing the initial position on 
an equipotential as r(0) = p e (</?o)(cos<^OjSin<po), we get 
t = ± J* o dip'J e {<p')/n. With 29, = f^dtpj t (<p) one 
obtains 

' M -t/r , M ^ TXW*(ip ,ipo±2n;e) 
dte t/T w^(tp(t);e) = it _ , (A4) 

where X — l/(rO) and 

W*(<p ,<pi',e)=T dipj e {ip)w ll {ip;e)e J *o 

(A5) 

To evaluate the average over initial values in Eq. (||), 
we first integrate along the equipotentials with fixed en- 
ergy e and then over e. It turns out that the Jaco- 
bian of the transformation from polar coordinates p, p 
to t he e nergy-angle coordinates e, tp is just given by 
Eq. (A3), dppdtp = de dtp J t (tp) . If the effective potential 
w(£, 77) is an even function of both arguments, we have 
Wp(ip + 7r;e) = -u^(<p;e) and J € (cp + n) = J e (ip), and 



all integrals over intervals of length 2ir can be reduced to 
integrals over intervals of length n, and we obtain 



xw p (v? ;e 



de / d<paJ e (ipo) 
Jo 

2XW^(p ,p ±7r;e) 

1 _l_ v-X8 e 



(A6) 



with the upper sign and a+ = e sup , 6+ = e max for 
equipotentials around a maximum and the lower sign and 
a_ = e m i n , 6_ = £i n f for those around a minimum. 

Open equipotentials of a periodic potential with rect- 
angular symmetry, by definition, connect one point on a 
boundary of the unit cell with the equivalent point on 
the opposite boundary. Since equipotential lines cannot 
cross each other, open equipotentials can exist either in x 
direction or in y direction, but not in both. Let us assume 
that in the energy interval ei n f < e < e S up open equipo- 
tentials in y direction exist. We may describe them in 
polar coordinates choosing the origin in a maximum, so 
that we can use the formalism developed above. Let the 
equipotential with energy e hit the upper boundary of the 
unit cell at r\ = ir = p e (tp e ) sintp e for 7r/4 < tp e < tt/2. 
Assuming that w(£,r/) is even with respect to both ar- 
guments, we can show that the GC motion in negative 
y direction with initial conditions tt — ip e < ipo < n + ip e 
yields the same contribution to the diffusion tensor as 
those moving in the positive y direction with initial con- 
ditions — tp t < tpo < ip e , and we consider here only the 
latter. 

For the time integration, we divide the infinite 
time interval into an initial one of duration t e = 
Q -1 J * dip J e ((p) and subsequent intervals of duration 

T e = SI -1 dp J e (ip), which is the time a GC needs 
to traverse a unit ce ll o n the equipotential of energy e. 
Using the definition (A5), we obtain for the contribution 
of open orbits to the diffusion tensor 

To this we have to add the contribution of closed orbits 



according to Eq. @, Df v = D%?* + D$g 
result can be written in the form of Eq. (|2 



D 



fMU 



The 



APPENDIX B: ANALYTIC AND ASYMPTOTIC 
RESULTS 

We present explicit results for the additive cosine 
model (20) of the effective potential, with < k < 1. 
The partial derivatives are then = — sin £ and = 
—k sin rj. 

We consider first the symmetric case n = 1 which, ac- 
cording to Eq. (||), can hold for all values of the magnetic 
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field only if the original modulation has square symme- 
try with equal lattice constants = 2n/K in 
both directions, so that Q = (2tt / a) 2 V x / (toq) • 

Things become especially simple close to the maximum 
at the origin, where w — £ and w n ~ —r\. Then the 
equipotentials become circles wit h ra dii p f = 4 — 2e inde- 
pendent of ip, and the Jacobian (A3) reduces to J~ e = 1. 
The angular velocity dip/ 1 dt — Q becomes constant along 
the equipotentials, and independent of e. Thus, the GC 
motion in this approximation is very similar to the simple 
cyclotron motion, however with the circular frequency Q 
instead of the cyclotron frequency ujq. As a consequence, 
all integrals in Eq. ( |A6| ) can easily be evaluated analyti- 
cally, with the result 



D cl ' + (e ) 



2 Q 2 T 7r(2-£ 



quj 



2tt7 (2tt) 2 1 + (Qt) 2 



(Bl) 



where e qu is the energy above which the quadratic ap- 
proximation is valid, and = D xx + , and D~£ + = 
-Df y + = QtD+ x . Thus, for Qt > 1, the motion of the 
GCs along closed equipotentials leads to a suppression 
oc (Slr)~ 2 . For a suitable choice of e qu (= 2-^/tt = 0.228) 
and a corresponding treatment of D^iT , one obtains the 
result indicated in Fig. || by the long-dashed line. 

Going beyond this simple quadratic approximation, we 
obtain qualitatively similar results. All equipotentials 
with energy e > (< 0) are closed lines around a maxi- 
mum (minimum). As |e| becomes small, the angular ve- 
locity varies along the orbits and becomes very small near 
the saddle points [(£,77) = (0,7r) and equivalent], where 
the Jacobian J £ {ip) diverges. Only the equipotentials ex- 
actly at e = are open trajectories (straight lines), but 
they yield vanishing contribution to the diffusion tensor. 
Exploiting the symmetry, we can show that D c ^ + and 
Dfy~ of Eqs. ( |A6D yield identical contributions to the di- 
agonal components of the diffusion tensor, whereas their 
contributions to the off-diagonal components cancel. The 
result for the non-vanishing diagonal components can be 
written as Eq. (U?). The numerically calculated $(f2r) 
is plotted in Fig. |. In the weak-modulation limit (or 
"dirty limit", fir — > 0) we have $(0) = 1, as expected 
from Eqs. ( |12| ) and (|15|). For Qt —> 00, $(Or) becomes 
small. We can expand Eq. (A6) for large Qt and show 
that the term linear in 1/Qt vanishes identically. The 
prefactor of the leading term can be calculated numeri- 
cally, and we obtain &(Qt) w 3.29/(Jlr) 2 for Sir — > 00. 
This can be used to obtain the one-parameter interpo- 
lation $ x (f2r) = 3.29/(3.29 + (fir) 2 ), which approxi- 
mates $(f2r) well for large values of Sir (see Fig. ||). 
An apparent improvement at small and intermediate Sir 
is obtained with the approximation $3 (Sir) defined in 
Eq. @. 

We now turn to the general rectangular symmetry. For 
< k < 1 the equipotentials with energies |e| < 1 — k are 
open (in y direction), and degenerate into straight lines 
for k = 0. 



In the "dirty limit" Sir — > the distinction between 
open and closed equipotentials is not relevant, since we 
can expand the velocity v^(t) into a Taylor series for 
small t and perform the integral in Eq. (|^) term by term. 
Up to second order in Sir we obtain 



$ xx (Qt,k) = k 2 [l-(Slr) 2 /2 + 

$yy(QT, K) = 1 - K 2 {Qt) 2 /2 + . . 



(B2) 
(B3) 



for all values of k. 

In the "clean limit" , open equipotentials dominate 
<& yy {Q,T, k) and introduce a characteristic \[k dependence 
for k 1 ■ Already the fraction of the unit cell covered 
by open equipotentials, ^ opcn /(27r) 2 = 1 - A^ e J(2w 2 ), 
which is plotted versus re (as dashed line) in Fig. y, shows 
such a dependence. To see that, we calculated the cor- 
responding area ^4™oscd = Si-k ^ e J-n & (</?) covered 
by closed equipotentials around a maximum (equal to 
that around a minimum), which allows the expansion 
^cios X cd = 16V^+0(re 3 / 2 ). 

The contributions of closed equipotentials to both <& xx 
and <&yy vanish in the clean limit. The contribution of 
open equipotentials to <j> TO (oo, re) is finite while that to 
^xx{oo,k) vanishes, because the average value of the 
guiding center velocity component v y (t) is finite, while 
that of v x (t) is zero. Since <& xx (Qt, k) behaves similar 
to & xx (Qt ) in the square-symmetric case, we extrapo- 
lated Eq. (P32| ) to arbitrary values of Qt and found that 
Eq. (E3) provides an extremely good approximation. 



<I> (k), clean, numerical 
O™(k)=1-1.65k' 2 +0.65k 5 ' 




FIG. 9. Numerical result in the clean limit for & yy (oo,K) 
(diamonds) versus anisotropy parameter k = V y /V x , to- 
gether with analytic approximation (solid line). Also shown 
is the fraction of the unit cell covered by open equipotentials 
(dashed). 

Since, in the clean limit, X = 1/(Qt) — > 0, we can 
easily evaluate Eq. flAE| ) along an open equipotential, 
W x (— ip e , ip e ; e) — —2ir, we obtain asymptotically from 
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Eq. @ *»»(oo,k) = de4/(fiT £ ) with 

T £ = 2K(g)/(^^), g - (4«/[(l + k) 2 - e 2 ]) 1 ' 2 , 

where K(g) is the complete elliptic integral!!. Ex- 
panding this for k 1, we obtain the leading terms 
& yy (oo,K) « 1 — (31/67t)k 3 / 2 . Adding a suitable term 
to satisfy $ TO (oo, 1) = 0, we obtained the approximation 



<& TO (oo, k) ss $2(«0 as defined by Eq. (|23|). Apparently 
the plot in Fig. ^ reveals only for n > 0.7 slight deviations 
between exact and interpolated result. 

Using Eq. ( |23| ) and the small fir expansion ([B3|), we 
tried to approximate <S> yy (yiT, k) by the one-parameter 
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